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We extend the numerical simulations of She et al. [Phys. Rev. Lett. 70, 3251 
(1993)] of highly turbulent flow with 15 < Taylor-Reynolds number Re\ < 200 up to 
Re\ ~ 45 000, employing a reduced wave vector set method (introduced earlier) to 
approximately solve the Navier-Stokes equation. First, also for these extremely high 
Reynolds numbers Re\, the energy spectra as well as the higher moments - when 
scaled by the spectral intensity at the wave number k p of peak dissipation - can be 
described by one universal function of k/k p for all Re\. Second, the ISR scaling 
exponents ( m of this universal function are in agreement with the 1941 Kolmogorov 
theory (the better, the large Re\ is), as is the Re\ dependence of k p . Only around k p 
viscous damping leads to slight energy pileup in the spectra, as in the experimental 
data (bottleneck phenomenon). 
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1 Introduction 



1.1 Universal turbulent spectra 

Very high Reynolds number turbulence still resists full numerical simulations. While 
in experiments Taylor- Reynolds numbers Re\ up to Re\ = 13 000 have been reported 
the most turbulent numerical flow has - to our knowledge - been realized by 
She et al. ||, who achieve Re\ = 200 at a resolution of 512 3 . These authors find 
that for all Re\ up to Re\ = 200 all their energy spectra coincide when scaled by 
the spectral intensity at the wave number k p of peak dissipation. I.e., the function 

(|u(k)|( 2 ))/(|u(k p )| 2 ) = F( 2 )(A ; A P ), (1) 

is universal as assumed by Kolmogorov and Obukhov - both in the inertial sub- 
range (ISR) and in the viscous subrange (VSR). Note that F^'^l) = 1 by definition. 
The universality is also found in experimental spectra 0. For further numerical 
simulations, see also ||. 

Kolmogorov and Obukhov || not only assumed universality of F^ 2 \k/k p ), but 
also its power law behavior F^ 2 '(k/k p ) oc (k/k p )~& in the ISR with the classical 
exponent £2 = 2/3. (The scaling exponent 2/3 in the discrete Fourier transform, 
which we use here, corresponds to 5/3 in the continuous case.) 

It has been argued in a long lasting debate (see e.g. |], |7|, §]) that there are small 
intermittency corrections to the classical scaling exponent (2 — 2/3. Unfortunately, 
even at today's state of the computational art M, ( 2 cannot sufficiently precisely 
be determined from full numerical simulations so that one could confirm or rule out 
deviations from 2/3. The reason is that for the tractable Re\ the available wave 
number range is quite narrow (k p /k min ~ 5 in 0; k min is the lowest wave number 
free of forcing). To obtain an ISR which extends over more than a decade, k p must 
be larger than 50/c m j n 0. To realize this, a resolution > 1500 3 is required 0. The 
required computer work increases as Re\ log 2 Re\ || . We are thus far away from 
being able to create developed turbulence in a numerical flow for as high Re\ as 
in experiment [Q. The huge gap between experiments and simulations is demon- 
strated in table 1. We therefore still need reasonable approximation techniques to 
numerically solve the Navier-Stokes equation. 

1.2 Reduced wave vector set approximation for high .Re- 
turbulence 

Such an approximation has been introduced by us in |JX0|, [TT|| . Meanwhile we could 
considerably improve our approach [|1^]. Here we employ it to study universal fea- 
tures of fully developed turbulence for Reynolds numbers Re between 730 and 1.4- 10 7 
(corresponding to Taylor- Reynolds numbers Re\ between 120 and 45 000, see table 
2 below). This leads to remarkably large ISRs. For Re — 1.4 ■ 10 7 , the extension of 
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the ISR is k p /k min « 2000, i.e., more than three decades, compared to k p /k min « 5 
achieved in the presently best Navier-Stokes simulation 0. 

For completeness we shortly repeat the main idea of our approximation. To deal 
feasibly with the many scales present in turbulent flow, we only admit a geometrically 
scaling subset K of wave vectors in the Fourier sum, K = Ujifj, thus Wj(x, t) = 
YskeK M i(k, t) exp (zk • x). K = {kffl, n = l,...,n max = 80} contains appropriately 
chosen wave vectors, K x = {k$ = n = l,...,n max = 80}, I = l,...,l max , 

are scaled replica of K Q . The choice of l max depends on the control parameter 
v, the viscosity. The incompressible Navier-Stokes equation is solved on K with 
periodic boundary conditions in a box of size (2irL) 3 . All lengths will from now 
on be measured in multiples of L, so the smallest component of the smallest wave 
vector is 1. The flow is permanently, non-stochastically forced on the outer length 
scale with energy input rate e. All times will henceforth be given in multiples of 
(L 2 /e) 1 / 3 , i.e., the energy input and thus in the stationary case also the energy 
dissipation rate is e = 1. The type of forcing does not influence our results sizeably 
||12j| . The smallest wave vectors whose amplitudes are free of forcing have length 
kmin — 3. The smallest wave vectors at all are ±(1, 1, —2) + permutations and thus 
have length a/6 = 0.82A; m j n . All wave vector amplitudes u(k, t) with |k| > k min = 3 
are free of forcing. 

1.3 Contents of the paper 

The paper is organized as follows: In section 2 we examine the Re- dependence of 
the Taylor-Reynolds number. In section 3 we confirm universality up to the highest 
Re we can treat. The form of the universal spectra is discussed in section 4. Section 
5 is left for a short summary of our findings. 



2 i?e-dependences 

We now come to our results. First, what are the Reynolds numbers Re of our 
approximate numerical turbulence? As usual, there is some arbitrariness in the 
definition of the Reynolds number. We regard L = A/2 as the outer length scale, 
where A is the wave length of the smallest wave vector, and Uq = 2u l rms as the 
typical velocity difference on the outer lengthscale. Thus Re = UqLq/v can be 
considered as an appropriate definition of the Reynolds number. The data for five 
simulations, covering four decades of Re, are given in table 2. We also list the Taylor 
Reynolds numbers Re x = ui irrns \ T ayior/v, where X Tay ior = ui,rms/{diUi) rms are the 
Taylor lengths. 

The dissipation rate is balanced by the input rate e ~ Uq/L |T3], We 
therefore write 

e = cJJl/L , (2) 
where c e is a dimensionless number. Since we choose e = 1, Lq = 7r/\/6, and since 
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we find Uq = 2-u ljrms from the numerical solution, we can determine c e from this 
equation. It turns out to decrease with increasing Re, seemingly to a final level 
somewhere near 6 • 10~ 3 . Note that for laminar flow, on the other hand, it holds 
c e oc Re -1 , see e.g. 

Equation (□) leads to the relation 



When c e eventually becomes universal, i.e., independent of Re, the well known 
large Re limit relation Re\ oc \f~Re~ is recovered. In experiment, for smaller Re the 
measured Taylor Reynolds number Re\ turns out to be smaller than predicted by 
Re\ oc \f~Re~ (15|, [Uj, as in fact c £ depends on Re. Both the large Re behaviour of 
Re\ with the power law exponent 1/2 and the deviations for smaller Re can also be 
seen in our approximate solutions, see Fig. [I]. For a more detailed discussion on the 



.Re-dependence of c e and Re\, see [17], [II 

One remark concerning the nominal value of the prefactor in eq. (^). Taking the 
large Reynolds numbers, we have (with c e = 6.5 ■ 10~ 3 from table 2) Re\ ~ 12\/Re, 
whereas experimentally it is Re\ ~ \fRe~, i.e., we overestimate the Taylor Reynolds 
numbers by one order of magnitude. We explain this as due to our approximation, 
as in our reduced wave vector set K the larger wave vectors are considerably thinned 
out. So less energy than ought to be is transported downscale, leading to a larger 
Ui^ rms than in real turbulence and thus to a larger Re\ = ui rms /i i {diUi) rrns . The 
same blocking effect by phase space sparseness at larger k also leads to an overestima- 
tion of the Kolmogorov constant b in the velocity structure function D(r) = 6(er) 2//3 



as already mentioned in |T2|, |T]J - also by a factor of about 10 (b m 70 instead of 



b = 8.4 as in experiment, see [0]), because both Re\ and D(r) are quadratic in 
Ultras- (Note that the denominator (diUi) rms is proportional to y/e and thus kept 
fixed.) 



3 Universal wave number spectra 

The wave number spectra of our approximate solutions are shown in Fig. |2|. We can 
confirm that F^(k/k p ) in fact is universal for all Re even much beyond Re\ = 200 
as studied in 0. As in [Q, the wave number k p of peak dissipation is found to 
increase as Re 3 ^, i.e., i]k p is constant also for the huge Re we simulated, see table 
3. We find k p ~ I/IO77. For Re = 1.4 • 10 7 the ISR extends as far as 3.3 decades to 
the left of k p , see table 2. In k-space \og 10 {k p /k min ) gives the extension of the ISR. 
This extension of the scaling range should also follow from the Reynolds number 
Re. In r-space the scaling occurs between the outer length scale Lq and 10rj, where 
rj = z/ 3 / 4 /e 1//4 is the Kolmogorov length scale |13|, [L4|]. For large Re it is 

L /10 V = c' v Re 3 /\ (4') 
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Formally eq. (4') can be derived from eq. (Q) and one gets = c^VlO. Note that 
thus c' v is also slightly -Re-dependent via c e (Re). Yet our c t is smaller than the 
experimental one, see last paragraph of section 2. Therefore we rather define a c v 
from the extension of the scaling range found in our numerical simulation than from 
(4'), namely by 

kmin/ kp C n Rc ^ . (4) 

This corresponds to (4') since L is approximately l/k min and IO77 is approximately 
l/k p . More precisely we have L = 3ir / (V6k min ) , so Lq 1 = 0.26k min . The ratio 
k p / k m i n can be extracted from our numerical spectra in k-space. For the largest Re = 
1.4- 10 7 we find kp/k min = 1950 and thus obtain c„ = (k p /k m i n )Re~ 3 ^ = 8.52 • 1CT 3 . 
We now disregard the small .Re-dependence of c v and calculate the extension of the 
scaling ranges for smaller Re from (^) with c v = 8.52 ■ 10~ 3 . The results are given 
in the last column of table 2 and are found in excellent agreement with the length 
of the scaling ranges seen in our approximate solutions for these smaller Re, see 
6th column of table 2. Only for the smallest Re we find disagreement, as expected, 
because for Re = 730 the constant c v in eq. (H) can no longer be considered as 
independent of Re. 

Moreover, also the F^(k/k p ) = (|u(k)| m )/(|u(k p )| m ), m=3,4,6,8,10, are found 
to be universal (see Fig. [| for m — 6). Due to the very extended ISR we can 
determine the power law exponents ( m of the universal functions F^ m \k/k p ) oc 
(k/k p )~^ m rather precise, cf. table 3. This is still not possible in full simulations, as 
the universal scaling range is too small, k v j k m i n < 5 0, Q . If a power law fit for full 
simulations is tried nevertheless, one gets scaling exponents ( m much smaller than 



the classical ones, due to the nonuniversal large scale forcing [19 



4 Form of the universal spectra 

Having shown universality, we now check several fits to the universal spectra 
F^ m \k/kp) that we obtained from our numerical data. 

4.1 Power law fit and intermittency 

Firstly, we fitted the spectra for all Re by the two parameter function 

F^ m \k) = k~^ exp [-(jfe - l)/k ( r ] }, k = k/k pj k^ = /k p . (5) 



This form has theoretical support [^TJ and was also successfully used to fit exper- 
imental spectra ^2|. The crossover between the power law behaviour in the 
ISR and exponential fall off in the VSR takes place at the wave number k^. In 
table 3 we list the parameters, obtained from a fit in the range < k < rj^ 1 /4 or 
< k/kp < 2.5, respectively. The fit (|5|) gives ( 3 rather near, but not exactly equal 
to 1 as it should be according to Kolmogorov's structure equation |14[]. This tiny 



deviation is corrected in table 3 by dividing the scaling exponents obtained from the 
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fit (Hi by £3. The resulting scaling exponents take their classical values ( m = m/3 
with high accuracy, see table 3. 

From the appearance of the classical scaling exponents ( m = m/3 one might 
deduce that there is no intermittency in our signal. This conclusion would not be 
correct. In fact, for small scales we do observe strong intermittency in the signal 



12fl . We therefore suggested to introduce local ( m (k), defined by local fits of 
type ([5D in the restricted k-ranges [k/y/W, k\^10] for all k, keeping the fixed at 



their global values [E|. 

The surprising result is shown in Fig. [| There are large intermittency correc- 
tions 5( m (k) = (m(k) — m/3 at small scales (VSR), only moderate intermittency 
corrections at large scales (stirring subrange, SSR), but hardly any deviations from 
classical scaling in the ISR. This astonishing result was extensively discussed already 
in [12|. Here it can be confirmed for a considerably larger .Re-range. 

In addition, we fitted eq. @ to our spectra, but now with k™ fixed at k^ = 
2k^ /m, fcjj = (13.5?7)~ 1 , see below. Again we find only tiny global intermittency 



corrections, which clearly decrease with increasing Re, as predicted by [15, E3f. For 



Re = 1.05 ■ 10 4 we find 5( 2 = 0.012, c^e = -0.058, going down to 5( 2 = 0.002, 
S(q = —0.011 for Re = 1.4- 10 7 , suggesting, that intermittency might be a finite size 



effect. For details and a theoretical explanation we refer to ref. [24 



We now discuss the crossover scale between ISR and VSR. From the fit (gD we 
get (k^)' 1 « 13.5m i.e., the crossover scale is one order of magnitude larger than 



the Kolmogorov length. This fact has long been known from experiment |TJj and 

theory p5| . From the maximum condition for k 2 F^(k) at k = k p it follows that 

~ (2) 

kj — 1/(2 — £2) = 3/4. In our simulations there are some fluctuations around 
this value. The reason is that we have discrete wave vectors which are not dense in 
the VSR. Thus k p can only be determined with limited accuracy, cf. Fig. [|, where 
the energy dissipation rate e{k) = k 2 (|u(k)| 2 ) is displayed. Typically the relative 
k-distances are 5k jk m 1/10, which corresponds to the ~ 10% deviations of 
from 0.75 in table 3. Of course, to increase the accuracy of k p , one could also define 
it in terms of from the global fit fl5|) with fixed (2 — 2/3 to be k p = Ak^ /3. The 
resulting small changes are not visible in Fig. 

Where does the crossover from ISR to VSR behaviour take place in higher order 
moments? From the fit (|5|) we find k^ = 2k^/m } i.e., k^ = 3/(2m) to a very 



high precision | I2| . This means that the ISR is considerably smaller for higher order 
moments, or, to state it differently, higher order moments feel viscosity earlier than 
lower order moments. 



4.2 Energy pileup in the crossover region between ISR and 
VSR 

We also applied fits different from as from Fig. |5] it might seem that (^) only 
badly fits the spectrum in the range around k p . The same observation was re- 
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ported already by She and Jackson || when they determined the universal function 
F^ 2 '(k/k p ) from experimental data. To improve the fit, they suggested to use the 
empirical three parameter function 

1-2/3 

F i2) (~k) = — (l + aW) exp [-{k - l)/kf% k = k/k p} kf = kf'/k p , 

(6) 

and found a = 0.8, (3 ~ 0.7, i.e., near k = k p the decay of the spectral power 
is diminished. Their physical interpretation is a pileup of excitation around k p , 
possibly due to coherent vortex structures. 

We also tried the fit (g) and found a 2, (3 « 1.8, and kjp' ~ 0.4, all slightly 
depending on Re. Thus the energy pileup at k « k p in our approximate Navier- 
Stokes solution seems to be even stronger than in experiment, as we have an additive 
correction term with a larger exponent. We ascribe the larger correction summand to 
our approximation: larger wave vectors are more and more sparse, so the dissipative 
scales cannot acquire the kinetic energy as fast as they should and pile it up around 
k p . Of course, F^(k/k p ) does not increase with k/k p as the correction term is 
strongly damped by exp (—k/k^') with kjp' much smaller as k^ before. The energy 
pileup can also be observed in Fig. [|, where we replotted F^ 2 '(k/k p ) with the two 
fits eq. ([]), where we fixed £ 2 = 2/3, and eq. (0). For k ~ k p the fit (^) is slightly 
superior to the fit (^). We have to raise doubts that the energy pileup is due to 
coherent vortex structures, as was speculated in [Q, because these are less in our 
approximation than in experiment, cf. [[K]], while our pileup is stronger than found 
in the measured data. 

Instead, the energy pileup may possibly be explained by the socalled bottleneck 



phenomenon |26]. This phenomenon can be described as follows: Imagine a triad 
Navier-Stokes interaction between the amplitudes u(ki), u(k2), and u(ks), ki + 
k 2 + kj = 0, k\ < &2,3, so that u(k 2 ) and u(ks) are already considerably damped 
by viscosity, in additon to the power law decrease oc k~ 2 ^ 3 . So the turbulent energy 
transfer downscale ~ kiu(ki)u(k2)u(k 3 ) would be reduced and stationarity could 
not be achieved, if u(ki) did not increase, i.e., an energy pileup at k\ is established. 
The effect is strongest if k\ is around k p , because there u(k 2 ) and u(k 3 ) are already 
considerably damped. Of course, there is also viscous damping around k\, which 
would counteract the bottleneck effect, but for k\ < rf x the damping by viscosity v 



is weaker than the damping by the eddy viscosity 1 2(| . 



For k <C k p Falkovich [^6j calculates the first order correction to the K41 spec- 
trum due to the bottleneck phenomenon and obtains 

F^Ck) oc 4 2 ii(^)[l + fc 4/ 7log(£ -1 )], k = k/kp. (7) 

We tried to compare our results for the scaling of the correction term with his 
prediction oc k 13 , (3 1 = 4/3 (apart from the log-correction), but as k k p (i.e., 
k <S 1) has to hold, we can only fit with a very limited number of data points. As 
in addition for small k the design matrix of the fitting problem nearly degenerates, 
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our results strongly depend on the fit range we choose and we can neither verify nor 
rule out his prediction (3' = 4/3. 

The same pileup as for F^ 2 \k) also appears in higher order velocity moments 
(|u(k)| m ) with m > 2, leading to smaller local scaling exponents ( m (k) for k near the 
VSR. Possibly this effect mimicks intermittency corrections to classical scaling in 
experimental data or simulated data with shorter scaling ranges which in fact would 
not show up in sufficiently long ISR for large enough Re. This interpretation is also 
supported by the behavior of the scale resolved intermittency exponents, Cm{k), see 
Fig. |3] above and ref. |TIJ . 



4.3 Log-similarity description of the spectra 

Finally, besides the normalization of the spectra (|I|) another procedure has been 
suggested to get an universal description of the experimental data, namely, the 
log-similarity description | 15fl . This claims that the logarithmic spectra coincide 



for different Re\, when both the abszissa and the ordinate are multiplied by some 
function /3(Re\), i.e., /31og(|u(k)| 2 /|u(k )| 2 ) against /31og(fc/fc ) is claimed to be 
universal. k is a wave number which has to be fitted to the experimental data. In 
0J fc « 2A; p « (677)- 1 is found. For large Re x it is (3 = 0.9/ log(i?e A /75) pj. We 



plot the spectra for our three largest Re in this parametrization, see inset of Fig. 
||. As Taylor Reynolds number we simply take \/Re, because our approximation 
overestimates Re\, see the discussion in section 2. For smaller Re the function 
P(Re\) behaves quite different, so it is not reasonable to show also the spectra for 
smaller Re in the plot. The quality of the superposition of the spectra might improve 
if one readjusted the free parameters of this description. 



5 Summary 

To summarize, we first demonstrated the universality of the normalized moments 
F^ m \k/k p ) up to Re = 1.4 • 10 7 within our approximate Navier-Stokes solution. 
Second, the F^ m \k/k p ) can pretty well be described by 

F^ m \~k) = k- m / 3 exp [~^(k - 1)), k = k/kp. (8) 

Near peak dissipation an additive correction to the classical scaling F^ m '(k/k p ) oc 
(k/k p )~ m / 3 shows up, which might be explained by an energy pileup around k v due to 
a viscosity induced bottleneck phenomenon p6f . This bottleneck phenomenon might 
mimick intermittency corrections to classical scaling in experimental or simulated 
data with less resolution and for smaller Reynolds numbers. 
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Tables 





Re 


Rex 


kp / k m in 


# of modes 


best experiments |jj 
best simulations 
REWA [here] 


1.7- 10 8 
sa 4- 10 4 
1.4- 10 7 


13 000 
200 
45 000 


« 10 000 

5 

« 2 000 


oo 
4- 10 8 
80 • 3 • 13 = 3120 



Table 1 

In the first two lines Re, Rex, the length of the scaling range k p /k m i n , and the 
number of contributing modes are compared for the most developed experimental 
and numerical turbulence, respectively. In the third line we give the same data for 
the largest Re of our reduced wave vector set approximation (REWA). 



V 


Imax 1 


Re 


Rex 


c e (Re) 




\og w (e v Re^) 


5 


10" 


-3 


5 




730 


122 


46 • 10" 3 


0.24 


0.08 


5 


10" 


-4 


7 


1.05 


10 4 


801 


15 • 10~ 3 


0.97 


0.95 


5 


10" 


-5 


9 


1.25 


10 5 


3590 


9.1 ■ 10" 3 


1.79 


1.75 


5 


10" 


-6 


11 


1.37 


10 6 


13 600 


7.0 ■ 10" 3 


2.48 


2.53 


5 


10" 


-7 


13 


1.40 


10 7 


44 800 


6.5 ■ 10~ 3 


3.29 


3.29 



Table 2 

Results from our approximate solutions of the Navier-Stokes equation for various 
v. Imax + 1 is the number of rescaled wave number replica K\. The definition of 
Re is Re = UqLq/u. Here L = tc/V6, and U = 2w l r . ms describes the velocity 
difference across the outer scale, being determined for each v from our numerical 
solution. (ui trms depends on v and increases from 1.42 to 2.73 for the v in the table.) 
Rex = ui trms \Tayior/v, as usual. The coefficient c e is calculated according to eq. @. 
In the last two columns the extension of the scaling range, log l0 (k p /k min ) (found 
from the numerical solution), is compared with that calculated from eq. (Q), see 
text. 
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V 




/■ 


[KpT]) 


\ K d V) 


~k [2) 


5 


10" 


-3 


f) 656 

V 7 . \j 'J \j 


1 997 


10.2 


1 1 8 


n 87 


5 


10" 


-4 


0.678 


1.986 


10.6 


13.6 


0.78 


5 


io- 


-5 


0.672 


1.975 


9.2 


13.7 


0.67 


5 


io- 


-6 


0.669 


1.994 


10.5 


13.3 


0.79 


5 


10" 


-7 


0.668 


1.990 


9.1 


13.3 


0.67 



Table 3 

The fit parameters to our approximate solutions of the Navier-Stokes equation 

(2) 

for the same v as in table 2. k p is the wave number with peak dissipation, k d the 
cut off from our fit fl5|), if we fit the spectrum in the interval [0, r/ _1 /4]. In the last 
column their ratio Vjp = k d 2 ^/k p is given. From the condition k 2 F^(k) maximal 

~ (2) 

at 1, the cut off should be k d = 3/4. 
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Figures 



Figure 1: Re\ (squares) and c e (circles) as functions of Re. Left ordinate is Re\, 
right ordinate is c e . For large Re the power law Re\ oc y/Re holds, for small Re 
there are deviations, as the dimensionless number c e in eq. (0) depends on Re. The 



dependence of c £ on Re strongly resembles the experimental one, cf. |15|, [18], [13] 



Figure 2: Universal spectra F^ m \k/k p ) for m = 2 (flatter) and m = 6 (steeper) 
for Re = 7.3 • 10 2 (triangles), Re = 1.05 • 10 4 (crosses), Re = 1.25 • 10 5 (squares), 
Re = 1.37 • 10 6 (pluses), Re = 1.40 • 10 7 (diamonds). The arrows indicate the 
smallest wave number free of forcing, k m i n , for the respective Re. The smallest wave 
number of allk e K is 0.82fc min , see section 1.2. The dashed arrow labels k m i n in 
the simulation by She et al. , the dotted arrow marks k m i n in the one of Vincent 
and Meneguzzi M. 



Figure 3: Scale resolved intermittency corrections —5( m (k) for m = 2,4,6,8,10, 
bottom to top. v = 5 • 10~ 7 , Re = 1.4 ■ 10 7 . The fit range is [k/y/W, k^/lO] for all k. 



Figure 4: Universal energy dissipation rate e(k)/e(k p ) versus k/k p for v — 5 • 10 ' 
i?e = 1.4 ■ 10 7 . Also shown are the fits resulting from (|5|) (dashed) and (^) (solid). 



Figure 5: Universal spectrum F^ 2 \k/k p ) for z/ = 5 • 10~ 7 , Re = 1.4 • 10 7 . The fits 
(H) (dashed) and (|D (solid) are compared, both with the same £ 2 = 2/3. Inset: 
Quality of universality in the log-similarity description for the second moments, 
see text. The symbols mean Re = 1.25 • 10 5 (squares), Re = 1.37 • 10 6 (pluses), 
Re = 1.40 • 10 7 (diamonds). On the abscissa we plotted (3 log w (k / (2k p )) , on the 
ordinate /31og 10 ((|u(A;)| 2 )/(|u(2fc p )| 2 )) with = 0.9/log(i2e A /75), cf. |] 
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